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Abstract 

We analyze xp dependence of the Drell-Yan transverse momentum broaden- 
ing in hadron-nucleus collisions. In terms of generalized factorization theorem, 
we show that the xp dependence of the transverse momentum broadening, 
A((7|,)(xi7'), can be calculated in perturbative QCD. We demonstrate that 
I^{qj){xp) is a good observable for studying the effects of initial-state multi- 
ple scattering and extracting quark-gluon correlation functions. 
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Parton multiple scattering is responsible for many interesting and important phenomena 
in high energy collisions involving nucleus, such as transverse momentum broadening, energy 
loss, as well as the nuclear suppression of quarkonium states. Understanding the mecha- 
nism of parton multiple scattering and its effects is crucial for making precise predictions 
of the nuclear suppression of quarkonium productions, which maybe a potential signal for 
the quark-gluon plasma in relativistic heavy ion collisions Parton multiple scattering 
can happen both at the initial state or at the final state. The Drell-Yan pair production in 
hadron-nucleus collisions provide an excellent place to study effects of initial state parton 
multiple scattering. The nuclear dependence in the Drell-Yan transverse momentum spec- 
trum for the large qt region has been studied in QCD perturbation theory [0,^. In terms 
of generalized factorization theorem in QCD 0], the effects of multiple scattering can be 
expressed in terms of multiparton correlation functions which are as fundamental as the 
parton distributions. However, these parton correlation functions are not well determined 
yet. 

In this letter, we derive the xp dependence of the Drell-Yan transverse momentum broad- 
ening. We show that A{q^){xF) can be used as a good observable to study the effects of 
initial-state multiple scattering and parton energy loss. In addition, it can be used as an 
excellent observable for extracting information on multi-parton correlation functions. 

Consider the Drell-Yan process in hadron-nucleus collisions, h{p') +A{p) — »• £~^£~{q) + X, 
where q is the four-momentum for the virtual photon 7* which decays into the lepton pair. 
p' is the momentum for the incoming beam hadron and p is the momentum per nucleon for 
the nucleus with the atomic number A. Let qT be the transverse momentum of the Drell-Yan 
pair, we define the averaged transverse momentum square as 

In Eq. (|l|), Q is the total invariant mass of the lepton pair with = q^. Since single 
hard scattering is localized in space, only multiple scattering (at least, double scattering) 
are sensitive to the nuclear size (or A^/^ type dependence). Therefore, in order to extract 
the effect due to multiple scattering, we introduce the nuclear enhancement of the Drell-Yan 

{Qt) as 

, (2) 

which is often called the transverse momentum broadening. The broadening of the trans- 
verse momentum square defined in Eq. @j should be sensitive to parton multiple scattering 
between nucleons inside a large nucleus. 

In a perturbatively calculable hard scattering process, having an extra scattering between 
physical partons is suppressed by a power of the hard scale . Therefore, multiple scattering 
in momentum space between physical partons correspond to an expansion in power series of 
l/Q^. In this letter, we limit ourselves to double scattering in momentum space. With only 
single and double scattering, the broadening of the Drell-Yan transverse momentum square, 
A(g|,), can be parameterized as 

A(g2) =^ + 5^1/3 ^ (3) 
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which is consistent to existing data P,^]. In Eq. bA^^^-term represents the contribution 
directly from the double scattering which is explicitly proportional to the nuclear size (oc 
A^^^), with A the atomic weight of the nucleus target. 

In Ref. 1^, Qiu and Sterman argued that the factorization theorem for hadron-hadron 
scattering should also be valid at the first non-leading power in momentum transfer, 
which is essential for systematically calculating the double scattering in QCD perturbation 
theory. According to this generalized factorization theorem, we can expand the numerator 
in Eq. (|^) as 



2 2 d'^hA 



dQ'^dq^ 



a,b 

^ a,b 

+4>a/A{x) 



+ 



(4) 



where ® represents convolutions over partonic momenta, and represents the terms that 
are suppressed by higher powers of l/Q^. In Eq. (^, C^'^\ C^'^\ and C''-^^ are perturbatively 
calculable hard parts. (pb/h{x') is the parton distribution of the beam hadron, and ipa/Aix) 
is the parton distribution in the nucleus normalized by the atomic number A. Ta/A{x) and 
Tb/h{x') are the four-parton correlation functions ||^ in the nucleus and the beam hadron, 
respectively. 

Because of the well-known EMC effect, as well as effects of nuclear shadowing and Fermi 
motion, the nuclear dependence of the 4>a/A is nontrivial. However, if we parameterize 
the effective nuclear parton distribution (pa/A into A" times corresponding nucleon parton 
distributions, we found a very small power of a for a wide range of x 0. Taking the 
EKS98 parameterization of nuclear parton distributions [1^ as an example, we found that 
a ~ ±(0.02 — 0.03) for the x-range relevant to this study, which is much smaller than 1/3 
for the A^^^-tjpe enhancements. Therefore, H^^ in Eq. (^, which is proportional to 4>a/A, 
should have a very weak nuclear dependence. Similarly, H^^^ in Eq. (^) also has a weak 
nuclear dependence. On the other hand, the nuclear parton correlation function Ta/A has 

an explicit dependence on the nuclear size (oc A^^^) 0|n[], so as the H^^^ in Eq. m. 

According to the factorization theorem §], the denominator in Eq. (|l|) can also be expand 
in terms of power series: 



d (i(3-^°^ 



a,b 



dQ' 




1 + 



where /i represents both renormalization and factorization scale. In Eqs. 
quantities are normalized by the atomic number A. Substituting Eqs. ( 
Eq. (|1]), we obtain 



) and (I 
and d) 



(5) 

I, all 
into 
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2\hA 



a 



(0) 



(6) 



In deriving Eq. (j^), we kept only terms up to 0{A^^^ /Q'^), and dropped all power correction 
terms without the A^^^ nuclear enhancement. Similarly, for a nucleon target, we have 



2\hN 
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(0) 

N 



On 



(7) 



Substituting above Eqs. and (|^ into our definition of the nuclear broadening of the 
transverse momentum square in Eq. (0), we derive 



H 



(0) 



H 



(0) 

N 



a 



(0) 



On 



+ 



^A 



(8) 



where we neglected terms of 0{A^ /Q"^). The first term in Eq. (P) should have a very weak 
nuclear dependence due to the fact that effective nuclear parton distributions have a very 
week A" dependence. Therefore, the first term contributes to the a-term in Eq. (|^). In 
addition, the first term should be numerically very small due to the cancellation between 
the two terms inside the bracket. On the other hand, the second term in Eq. @ represents 
the double scattering contribution and it gives the A^/^-type enhancement, and therefore, it 
contributes to the hA^/^-teim in Eq. (^. 

At the leading order, the inclusive cross section is 



^A 



da 



hA^e+e- 



^oJ2'^qJ dx' (f)g/hix') J dx(t)qjA{x)5{Q'^ - xx's) 



(9) 



with s = {p + p') and the Born cross section 



0^0 



Also at the leading order, the Drell-Yan differential cross section is given by 

d<^hA^i+l~ d<^hA^l+l- 



dQ'^dqj' 



dQ^ 



6{q. 



T) 



(10) 



(11) 



Therefore, H^^^ and ifjj''' in Eq. 



vanish at the leading order, because / dq^q^6{q^) = 0, 
and consequently, the first term in Eq. (H), or the a-term in Eq. @, vanishes at the leading 
order in a^. 

The double scattering contribution to the Drell-Yan transverse momentum broadening, 



H^^'^ in Eq. 



H 



(2) 



was derived in Ref. |r3|. At the leading order in a^, it is given by \T3 



^ J dx'^g/h{x') J dxT^)\{x)5[Q^ 



XX s] 



(12) 



where the four-parton correlation function is defined as 
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X I (PA|F.+ (l/2-)^.(0)7"'^,(2/-)i^+"(yr)bA) . (13) 

These parton correlation functions are not well measured yet. By comparing the operator 
definitions of the correlation functions and the definitions of the normal twist-2 parton 
distributions, Luo, Qiu, and Sterman (LQS) proposed the following model |5|,pJ]l: 



Tf/A{x) = X'A'/'4>f/A{x) , (14) 



where A is a free parameter to be fixed by experimental data, and was estimated in Ref. |TB 
as = 0.01 GeV^ by using the Drell-Yan data from NAlO and E772 experiment 



In order to obtain the leading order x^-dependence of A{q^), we multiply 6{xf — {x' — 
x)) dxp to the right-hand side of Eq. (p!2D, and obtain 



dHf{xF) ( 47r^a, ^ ^ (/) ^ 1 



with 



Xi = (a/x|, + 4r + Xi?)/2 , and X2 = (-\/x|, + 4r — Xi?)/2 , (16) 



where r = Q'^/s. Combining Eqs. (|), (P), and (|l^), we obtain the xp dependence of the 
Drell-Yan transverse momentum broadening A{q^) at the leading order in 

dA{q^j,){xp) ^ fAfoA . e^0g-/fe(xi)T^^;i(x2)/(a:i + X2) 
dxp V 3 / Eg elJdx'(f)q/h{x')(f)g/A{r/x)/x' 

From Eq. ([17|), we see that dA{q'^) /dxp directly depends on the quark-gluon correlation 
functions. Therefore, the broadening is a very good observable for measuring quark-gluon 
correlation functions. 

If we use LQS model for Tg/A given in Eq. (|T^, from Eq. (0), we can derive a much 
simpler expression for the broadening, 

d^Qr) ^ ( X 2 .1/3 da j da_ 
dxp V 3 ; dQ^dxp/ dQ^ ' ^ ' 

From Eq. ([I8|), we can see that the Drell-Yan transverse momentum broadening should have 
similar xp dependence as the differential cross section da/dQ^dxp if the LQS model for 
Tg/j^(x) is valid. Therefore, by comparing the xp dependence of dA{q^) / dxp and da/dQ'^dxp, 
we can provide an immediate test of LQS model for the correlation functions. We empha- 
size that even if LQS model is not a good approximation for the quark-gluon correlation 
functions, measuring the x^-dependence of the nuclear broadening provides excellent infor- 
mation for extracting the quark-gluon correlation functions Tq/A{x) directly, as shown in 
Eq. (0). 

In the following, we use Eq. (|18D to obtain the numerical estimates of the xp depen- 
dence of the Drell-Yan transverse momentum broadening. Although the value of A^ for the 
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correlation function is not well determined, a different value of corresponds to a simple 
adjustment to the overall normalization. Therefore, the uncertainty in the value of should 
not affect our following discussions. 

In obtaining our following numerical results, we use the CTEQ4L distribution as the 
quark distributions in the nucleon. For effective quark distributions in the nucleus, we 
define qi/A{,x) = qi/p{x)Ri{x, A), and use EKS98 for the parameterizations of Ri{x,A) [jlO|, 
which fit the data well. We choose the renormalization and factorization scale to be /i = Q, 
and choose the incoming beam energy p = 800 GeV which is the energy used by the Fermilab 
experiments [J7tp!3|-p7|] . 

In Fig. 0, we plotted dA{q^) /dxp in Eq. (p^) as a function of xp with A = 184 and Q = 
5 GeV and 11 GeV, respectively. The dotted lines correspond to EKS98 parameterizations 
of effective nuclear parton distributions. In order to separate the nuclear dependence caused 
by multiple scattering and that caused by the effective nuclear parton distributions, we also 
plotted dA{q^) /dxp in solid curves with Ri{x, A) = 1. The difference between the solid and 
dotted lines is a direct consequence of the difference between nucleon and nuclear parton 
distributions. As shown in Eq. (p^Tj), X2 represents the momentum fraction of a quark (or 
antiquark) from the nuclear target. At xp = 0, we have X2 = Ql\fs ~ 0.13 for Q = 5 GeV, 
and ^ 0.28 for Q = W GeV. It is clear that Q = 5 GeV and Q = W GeV cover very 
different range of x^- For Q = 5 GeV, xp > ^ corresponds to 0:2 < 0.1 or corresponds to 
the shadowing region, while xp <^ covers the region of EMC effect. On the other hand, 
the Drell-Yan pairs of Q = 11 GeV are not sensitive to the nuclear shadowing at all. The 
entire range oixp values matches the range of EMC effect which include the antishadowing 
region (0.1 < X2 < 0.2), and the EMC suppression region (0.2 < 0:2 < 0.7), as well as the 
Fermi motion region for larger x^- Such dependence in effective nuclear quark distributions 
are clearly shown in Fig. In Fig. |I]a, the dotted line is above the solid line in the central 
region oixp due to the fact that X2 is in the antishadowing region. When Xi;' is positive, the 
dotted line is below the solid line because qijA^^'i) is now in the shadowing region. When xp 
is less than zero, the dotted line is again below the solid line due to the fact that X2 is now in 
the region of EMC suppression. On the other hand. Fig. |l|b shows slightly different relation 
between the dotted and solid line due to the fact that at Q = 11 GeV, X2 covers a different 
range as shown in Fig. [l|. Although the nuclear dependence in effective nuclear parton 
distributions change the xp dependence of the Drell-Yan transverse momentum broadening, 
the change (e.g., the difference between the dotted and solid lines in Fig. |l|) is extremely 
small. Therefore, the xp dependence of the nuclear broadening, A{q^){xp) is not very 
sensitive to nuclear shadowing, and can be an excellent probe of parton multiple scattering. 

In addition to the transverse momentum broadening, parton multiple scattering also 
causes the energy loss. In recent experiments, the ratio of the cross section per nucleon from 
the Drell-Yan pair production m p — A collisions has been used as a direct measurement 



of the parton energy loss in nuclear medium However, as pointed out in Ref. |T^, a 



substantial fraction of the variation in the cross section ratios versus Xi comes from the 
shadowing of 4'f/A{x2) at small X2, and therefore it is difficult to extract precise information 
on parton energy loss from the cross section ratios. On the other hand, BDMPS showed that 



for a jet produced in nuclear matter, the energy loss per unit length — ^ and the transverse 



momentum broadening /\pj_ has the following relation JTS 



6 



dE 

dz 



(19) 



Similar to Eq. ([T9|), one can use the transverse momentum broadening of the Drell-Yan pair 
A{q'^){xF) to estimate the parton energy loss due to the initial state multiple scattering. As 
we demonstrated in Fig. [l|, the variation in the shape of A{q'^){xF) due to the shadowing is 
small. Therefore, A{q'^){xF) is also a good probe for the parton energy loss. 

Cross sections on nuclear targets are often parameterized as A" times corresponding 
cross sections on a nucleon target. In principle, the power a can be a function of A, qt, 
xpi as well as other physical observables. The Nuclear dependence of a is often used to 
study the effects of multiple scattering. In recent experimental studies of the J/\l/ and 
suppression m. p — A collisions , data was presented in terms of a as a function of the 



transverse momentum px in different xp regions. It was found |TH| that the shapes of a{pT) 



are very similar in different regions of x^, and it was concluded |jT6[ that the parton energy 
loss is independent of the cc energy. 

Using our result of A{q'^){xp), we can also estimate the a as a function of qp in different 
xp region for the Drell-Yan production. We define the parameter a{qT) for the Drell-Yan 
production as 



dQ'^dqj. 



da 



hN 



dQ'^dqj- 



(20) 



If qx is not too large, the qp spectrum of the Drell-Yan pairs can be approximately param- 
eterized by a Gaussian form [|T^] , 



da 



hN 



dQ'^dqj, 



oc 



2^hN 



IT 



Qp 
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1 dahA 

A dQ^dq^ '^^ {ql) 
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2\hA 



exp 
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qp 



(21) 



2 \hA 



(22) 



Substituting Eqs. (0) and (p2D into Eq. (^OD, we derive 



a(gi 



1 + 



\n{A) 



In 



X 



qp 



2\hN 



(23) 



where x = ^{Qp) / (Qp)^'^ ■ In deriving Eq. (^), we used Eq. (^. In order to estimate the 
value of a{qp), we use the value (qp)^^ = 1.2 GeV^ for the cross section per nucleon [ll7| . 
And we use Eq. (|T8|) to integrate over different ranges of Xi? value to obtain the A{q^) in 
different ranges of xp. For A = 184, we choose three different xp ranges, which are the 
same as the regions used in Ref. [|l^], small xp (SXF: —0.1 < xp < 0.3), intermediate xp 
(IXF: 0.2 < Xp < 0.6), and large xp (LXF: 0.3 < xp < 0.93). For these three xp regions, 
we obtain the corresponding values of A{q^) as [0.10, 0.043, 0.026] GeV^ for Q = 5 GeV, 
and [0.078, 0.043, 0.027] GeV^ for Q = 11 GeV, respectively. It is clear that the value of x 
should be very small, and the power a{qp) in Eq. (|2^) can be approximated as 
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From Eq. (pi]), we conclude that if qt is not too large, a^qx) should show linear depen- 
dence on g|. (or quadratic dependence on g-j-). For the Drell-Yan production, we expect an 
extremely small coefficient due to a weak transverse momentum broadening (or a small x 
value in Eq. (||)). 

In the above discussion, we used the same averaged value of {q^)^^ for different xp 
ranges. In principle {q^)^^ should also have xp dependence due to kinematics. The larger 
\xf\, the smaller (q^)'^^- Hence {q^)^^ should have similar xp dependence as ^{q^){xp) 
shown in Fig. [1|. Therefore, x = ^{Qt) / {Qt)^^ should have even smaller dependence on xp. 

In Fig. H, we plotted our predictions for the value of ai^qp) defined in Eq. (p3|) as a 
function of qp in three different xp regions. From Fig. 0, it is clear that the shape of 
a{qT) in different xp ranges are similar to what was observed in recent data on J/\l/ and 
suppression [|TB|. Such similarity is natural because both qp spectrum of the Drell-Yan and 



Charmonium production can be approximated by a Gaussian form when qp is not too large. 
In addition, we emphasize that as shown in Figs. |T] and ^, even though a{qp) seems not to 
be sensitive to xp, the xp dependence of transverse momentum broadening, A{q^){xp), can 
be very strong. 

In summary, we derived the xp dependence of the Drell-Yan transverse momentum broad- 
ening in terms of quark-gluon correlation functions. We predicted that A{q'^){xp) has a 
strong Xp dependence. In particular, if the LQS model for the correlation functions is valid, 
A{q^){xp) should have similar xp dependence as the differential cross section da / dQ'^dxp. 
We also demonstrated that A{q^){xp) is a very sensitive observable to study the effects of 
initial state multiple scattering and the parton energy loss. In fact, A{q^){xp) itself is an 
excellent observable for extracting direct information on the nuclear quark-gluon correlation 
functions. 
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FIG. 1. The transverse momentum broadening of the Drcll-Yan pair, dA{q^) /dxp, as a function 
of xp for 800 GeV proton beam on nuclear target A=184, at Q = 5 GeV (a) and Q = 11 GeV (b). 
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Fig.2 

FIG. 2. The ai^qx) as a function of qt for smah, intermediate, and large xp region, respectively. 
The curves are for Q = 5 GeV and 800 GeV beam energy. 
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